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The field equations of a special class of tetrad theory of gravitation have been applied to 
t-h . tetrad space having three unknown functions of radial coordinate. The spherically symmetric 
^ .vacuum stress-energy momentum tensor with one assumption concerning its specific form 
7— \ [ generates two non-trivial different exact analytic solutions for these field equations. For 
^ | large r, the exact analytic solutions coincide with the Schwarzschild solution, while for small 
O r, they behave in a manner similar to the de Sitter solution and describe a spherically 
^ 'Symmetric black hole singularity free everywhere. The solutions obtained give rise to two 
different tetrad structures, but having the same metric, i.e., a static spherically symmetric 
q-i nonsingular black hole metric. We then, calculated the energy associated with these two 
^ | exact analytic solutions using the superpotential method. We find that unless the time- 

space components of the tetrad go to zero faster than —= at infinity, the two solutions give 
> . y/r 

; different results. This fact implies that the time-space components of the tetrad must vanish 
H ! faster than —= at infinity. 
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1. Introduction 



In the tetradic theories of gravitation the gravitational field is described by a tetrad 
field. In the past, several authors have considered these theories to solve some problems of 
general relativity or to extend it to include other fields. While searching for a unification 
of gravitation and electromagnetism Einstein considered a tetrad theory 0, but soon he 
rejected it because there was no Schwarzschild solution in his field equation. 

M0ller has shown that it is impossible to find a satisfactory expression for the energy- 
momentum complex in the framework of Riemannian geometry p|. In a series of papers, 
HI H M0ller was able to obtain a general expression for a satisfactory energy-momentum 
complex in the absolute parallelism space. The Lagrangian formulation of this tetrad theory 
of gravitation was first given by Pellegrini and Plebanski ||. In these attempts the admissible 
Lagrangians were limited by the assumption that the equations determining the metric 
tensor should coincide with the Einstein equation. M0ller || abandoned this assumption 
and suggested to look for a wider class of Lagrangians, allowing for possible deviation from 
the Einstein equation in the case of strong gravitational fields. Sdez generalized M0ller 
theory into a scalar tetrad theory of gravitation. Meyer || showed that M0ller theory is a 
special case of Poincare gauge theory 0, |T0 . 



Hayashi and Nakano JTT| formulated the tetrad theory of gravitation as the gauge theory 
of space-time translation group. Hayashi and Shirafuji [0 studied the geometrical and ob- 
servational basis of the tetrad theory, assuming that the Lagrangian be given by a quadratic 
form of torsion tensor assuming the Lagrangian to be invariant under parity operation in- 
volving three unknown parameters to be fixed by experiments, besides a cosmological term. 
Two of the three parameters were determined by comparing with solar-system experiments, 



1^1 while only an upper bound has been estimated for the third one [H 

The numerical values of the two parameters found were very small consistent with being 
equal to zero. If these two parameters are equal to zero exactly, the theory reduces to the 
one proposed by Hayashi and Nakano JTT] and M0ller, || which we shall here refer to as 
the HNM theory for short. This theory differs from general relativity only when the torsion 
tensor has nonvanishing axial- vector part. It was also shown pi that the Birkhoff theorem 



can be extended to the HNM theory. Namely, for spherically symmetric case in vacuum, 
which is not necessarily time independent, the axial- vector part of the torsion tensor should 
vanish due to the antisymmetric part of the field equation, and therefore, with the help 



of the Birkhoff theorem [|1J] of general relativity we see that the spacetime metric is the 
Schwarzschild. 



Mikhail et al. [15] derived the superpotential of the energy-momentum complex in the 
HNM theory and applied it to two spherically symmetric solutions. It was found that 
in one of the two solutions the gravitational mass does not coincide with the calculated 
energy. This result was extended to a wider class of solutions with spherical symmetry 



rj]. An explicit expression was given for all the stationary, asymptotically flat solutions 
with spherical symmetry, which were then classified according to the asymptotic behavior 
of the components of (A o) and (A a ). It was found that the equality of the gravitational and 
inertial masses holds only when (A o ) and (A a ) tend to zero faster than 1 / y/r. 
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Dymnikova [ljj derived a static spherically symmetric nonsingular black hole solution in 
orthodox general relativity assuming a specific form of the stress-energy momentum tensor. 
This solution practically coincides with the Schwarzschild solution for large r, for small r 
behaves like the de Sitter solution and describe a spherically symmetric black hole singularity 



free everywhere [17 



The general form of the tetrad, A M , having spherical symmetry was given by Robertson 

i 

lSfl. In the Cartesian form it can be written as R 



A° = %A, A° = Cx a , X a = iDx a 
o a o 

\ a = S^B + Fx a x a + e aoiP Sx p , (1) 

where A, C, D, B, F, and S are functions of t and r = (x a x a ) 1 ^ 2 , and the zeroth vector A M 

has the factor i = y/—l to preserve Lorentz signature. We consider an asymptotically flat 
space-time in this paper, and impose the boundary condition that for r — > oo the tetrad (1) 

approaches the tetrad of Minkowski space-time, (^^j — diag(i, S a a ). Robertson has shown 

that: 

1- Improper rotation are admitted if and only if S = 0. 

2- The functions C and F can be eliminated by a mere coordinate transformations, i.e., by 
making use of freedom to redefine t and r, leaving the tetrad (1) having three unknown 
functions in the Cartesian coordinate, which will be used in section 3 for calculations of 
HNM field equations but in the spherical polar coordinate. 

It is the aim of the present work to find the, asymptotically flat solutions with spherical 
symmetry which is different from the Schwarzschild solution in the HNM theory. Assuming 
the same form of the stress-energy momentum tensor as given by Dymnikova |T7| , we obtain 
two different exact analytic solutions of the HNM theory. We then, calculate the energy of 
those solutions using the superpotential of Mikhail et al. jjTJJ. 

In section 2 we briefly review the tetrad theory of gravitation. In section 3 we first 
apply the tetrad (1) with three unknown functions of the radial coordinate in spherical polar 
coordinates to the field equations of HNM theory, then we derived the corresponding partial 
differential equations. Assuming some conditions on these partial differential equations we 
obtained two different exact asymptotically flat solutions with spherical symmetry. In section 
4 the energy of the gravitating source is calculated by the superpotential method. The final 
section is devoted to the main results and discussion. 

Computer algebra system Maple V Release 4 is used in some calculations. 



*In this paper Latin indices ■■■) represent the vector number, and Greek indices (fx, v, ...) represent 
the vector components. All indices run from to 3. The spatial part of Latin indices are denoted by (a, b, ...), 
while that of greek indices by (a, (3, ...). In the present convention, latin indices are never raised. The tetrad 
A M is related to the parallel vector fields of @ by A^ 1 = ib^ and A M = b a ^. 
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2. The tetrad theory of gravitation 



In this paper we follow Miller's construction of the tetrad theory of gravitation based on 
the Weitzenbock space-time. In this theory the field variables are the 16 tetrad components 
A M , from which the metric is derived by 

i 

i i 

The Lagrangian L is an invariant constructed from 7 W and g^, where 7 W is the contorsion 
tensor given by 

l P Vp — A/jAi/;p, (3) 

where the semicolon denotes covariant differentiation with respect to Christoffel symbols. 
The most general Lagrangian density invariant under parity operation is given by the form 

C d ^ {-g) 1 ' 2 + a 2 'f up j llup + a 3 Y up l P v») , (4) 

where 

g d ^ f - det(^), (5) 
and $ M is the basic vector field defined by 

*, = W (6) 

Here ot\, a 2 , and a 3 are constants determined by M0ller such that the theory coincides with 
general relativity in the weak fields: 

ai = --, a 2 = -, ct 3 = -(l-A), (7) 

K K K 

where k is the Einstein constant and A is a free dimensionless parameter^. The same choice 
of the parameters was also obtained by Hayashi and Nakano |[TT| . 



M0ller applied the action principle to the Lagrangian density (4) and obtained the field 
equation in the form 

G pv + H P v = —kT^v, (8) 
F pv = 0, (9) 
where the Einstein tensor G^ v is defined by 

G P v = Rpu — -j9iivR- (10) 



'Throughout this paper we use the relativistic units, c = G = 1 and k = 8ir. 
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Here H„ u and F„ u are given by 



n def. , 



'ir. 



and 



®H,v ~ $»,p 



/ii/ i vp i • w -p 

and they are symmetric and skew symmetric tensors, respectively. 



(12) 



M0ller assumed that the energy-momentum tensor of matter fields is symmetric. In 
the Hayashi-Nakano theory, however, the energy- momentum tensor of spin- 1/2 fundamental 
particles has nonvanishing antisymmetric part arising from the effects due to intrinsic spin, 
and the right-hand side of (9) does not vanish when we take into account the possible effects 
of intrinsic spin. Nevertheless, since in this paper we consider only solutions in which the 
left hand side of (9) is identically vanishing, we refer the tetrad theory of gravitation based 
on the choice of the parameters, (7), as the Hayashi-Nakano-M0ller (HNM) theory for short. 

It can be shown 



12| that the tensors, and F^ u , consist of only those terms which are 



linear or quadratic in the axial- vector part of the torsion tensor, a„, defined by 



def. 



-e 



pupcr 



1 



(13) 



where en UpcT is defined by 



def. 



-p,Vp(7 



(-9) 1/2 *, 



(14) 



with 5n Up(T being completely antisymmetric and normalized as #0123 



1. Therefore, both 
is found to vanish 

from the antisymmetric part of the field equations, (9), the symmetric part (8) coincides 
with the Einstein equation. 



J pupa 

H^y and vanish if the a M is vanishing. In other words, when the a M 



3. Spherically symmetric nonsingular black hole solutions 



The tetrad space having three unknown functions of radial coordinate with spherical 



symmetry in spherical polar coordinates, can be written as [18 



A" 

i 



( iA iDr 

B sin 9 cos ( 

B sin 9 sin < 

Bcos9 





— cos 9 cos ( 
r 

— cos 9 sin < 

r 

sin 9 
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\ 

B sin (p 
r sin 9 

B coscj) 
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(15) 



5 



I. 



where % = 

Applying (15) to the field equations (8), (9) we note that the two tensors and 
are vanishing identically regardless of the values of the functions A, B and D. Thus M0ller 
field equations reduce for the tetrad (15) to Einstein's equations. Then the field equations 
(8) and (9) take the form 



00 



rA 2 B A 



3D 2 + 8B' 2 



^jD-2 (2DB" + B'D^j B^r 3 B 2 D 



2^DB" + B'D^jB- 5DB' 2 ^r 5 D 3 - (2BB" - 3D 2 - 3B ,2 ^jrB 4 + 



01 



2 \BD' - ADB'j r A BD 3 + 2 \BD' - QDB'j r 2 B 3 D - AB 5 B' 

D 

AB 1 



I^^DB" + B'D'^JB - 5DB' 2 ^r 3 D + ^2BB" - 3D 2 - 3B' 2 ^jrB 2 - 



2 1 BD' - ADB'j r 2 BD + AB 3 B 
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rAB A 



^3D 2 + B' 2 ^j A + 2BA'B'^rB 2 - ^2 ^DB" + B'D'^j B - 5DB ,2 ^r 3 AD + 



22 



2 \BD' - ADB'J r 2 ABD - 2AB 3 B' - 2B 4 A' 

J (dA" + 3A'D^j B - 3DA'B'^ABD + | ^2DB" + 5B'D'^j BD 



A 2 B 4 

( DD" + D' 2 \B 2 - 5D 2 B' 2 \A 2 - 2B 2 D 2 A' 2 \r 3 



B' 2 - 3D' \A Z - AB' Z A" - B"BA l + 2B Z A' 2 WB 2 - 



12 A 2 



2 a" d"D/|2 1 od2/|/2I„d2 



2| ^3BD' - ADB'^jA - 2BDa\t 2 ABD + A 2 B 3 B' + AB 4 A' 

T 33 = sm9 2 T 22 , (16) 

.. dA dB dD 

where A = — — , B = — — and D = ——. 

dr dr dr 

Now we are going to find some special solutions to the partial differential equations (16), 
assuming that the stress-energy momentum tensor has the form JET] 



T ° = Tl 1 , 
T 2 = T, 3 . 



(17) 



A first non-trivial solution can be obtained by taking D(r) = 0, and solving for A(r) and 
B(r), then we obtain 
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.4 



B 



1 



'i-^(i- e - fi v; 



where is a new radial coordinate defined by R = r / B and 



= r g r 2 , 
r g = 2m, 



r 



The form of the energy-momentum tensor is 

To = Tx 1 = coe-*/**, 
T 2 2 = T 3 3 = e e- fl3 /- 3 



3it! 3 \ 
2n 3 J ' 



(19) 



(20) 



and the tetrad (15) takes the form 



A' 



sin 9 cos (j)J 1 — (l — e 



-R 3 / ri 3 



R 



i 2m / 
sin6'sin0Wl — ( 1 — e 



-R 3 /ri 3 



R 



cos^l-^(l-e-^ 3 



with the associated Riemannian metric 



cos 9 cos < 
R 

cos 9 sin ( 
R 

sin 9 



R 



sin ■ 



Rsm9 



cos <p 
Rsin9 



(21) 



dR 2 



d s A = - Vl df + + R 2 dQ 2 
Vi 



(22) 



where 



and dVt 2 = d9 2 + sin 2 
solution W71 



2m 
~R 



(23) 



2 , which is a static spherically symmetric nonsingular black hole 
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A second non-trivial solution can be obtained by taking A(r) = 1, B(r) = 1, D(r) ^ 
and solving for D(r). In this case the resulting field equations of (16) can be integrated 
directly to give 



D(r) 



1 2m 
73" 



J g— r 3 /ri a 



Substituting for the value of -D(r) as given by (24) into (15), we get 



-r 3 /ri' c 









sin 9 cos (f) 
sin 6 sin 
cos6> 



cos 6* cos sin . 



r 

cos 6 sin < 

r 

sin 6 



r sin 

COS0 

r sin 




(24) 



(25) 



with the associated metric 



ds 2 = - 



1- — (l-e-^ 13 ) 



/2m 



■ x _ e - r y ri ^ drdt + dr 2 + r ^n 2 , (26) 



it is to be noted that m in the metric (22) and (26) is a constant of integration that will play 
the role of the mass producing the field in the calculations of the energy. Also the form of 
the stress-energy momentum tensor for this solution is the same as given by (20). 

Using the coordinate transformation 



Dr 

dT = dt + T^dr, 

we can eliminate the cross term of (26) to obtain 

dr 2 



(27) 



ds z = ~7] 2 dT z + + r W, 

V2 



(28) 



where r] 2 is defined by (23) and R = r. 

Thus we have two exact solutions of HNM field equations, each of which leads to the 
same metric, a static spherically symmetric nonsingular black hole in the spherical polar 
coordinate. 

The solutions (18) and (24) are the exact solutions of the HNM field equations. They 
practically coincide with the Schwarzschild solution for r >> rl and behave like the de Sitter 
solution, for r << rl. 

As is clear from (20), the spherically symmetric stress-energy momentum tensor is really 
anisotropic. The difference between the principle pressures 



n k = - Pk , 



(29) 
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correspond to the well known anisotropic character of evolution of the space-time inside 
a black hole undergoing a spherically symmetric gravitational collapse For r << t\ 

isotropization occurs and the stress-energy momentum tensor takes the isotropic form 



T a p = eg a p. (30) 

When r — > the energy density tends to e . For r >> r\ all the components of the stress- 
energy momentum tensor tend to zero very rapidly. 

The important result obtained in this section is that we have been able to derive two dif- 
ferent solutions for HNM theory; the Riemannian metric associated with these two solutions 
are identical, namely spherically symmetric nonsingular black hole. Since HNM theory is a 
pure gravitational theory, the above two solutions have to be equivalent in the sense that they 
describe the same physical situation. In what follows we examine the equivalence of these 
solutions by calculating the energy associated with each of them, using the superpotential 



derived for Miller's theory by Mikhail et al. [15 



4. The Energy Associated with each Solution 



The superpotential of the HNM theory is given by Mikhail et al. [y| as 



U» x = K -^—P XP a TuX [^g ax g,r - \g T ,i xp(T - (i - 2\)g T ^ px ] , (31) 



where P X pa TI/ is 



tvX del. r-r vX i r T „ vX X t„ vX 



p XP o = ^9 P r + s P T g« x v * - 5/g x r (32) 

with g pa uX being a tensor defined by 



9pa »* d ^ 5 p "6 a x - VV- ( 33 ) 



The energy is expressed by the surface integral [20 



E = lim / U 0a n a dS, (34) 

r ~ > °° J r=constant 

where n a is the unit 3-vector normal to the surface element dS. 

Now we are in a position to calculate the energy associated with the two solutions (18) 
and (24) using the superpotential (31). As is clear from (34), the only components which 
contributes to the energy is Uo° a . Thus substituting from the first solution (18) into (31) we 
obtain the following non-vanishing value 

W 0Q = ^-(l-e-^ 13 ). (35) 
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Substituting from (35) into (34) we get 

E{r) = m (l - e - r ' A/rl3 ) . (36) 

This is a satisfactory result and should be expected. 

Now let us turn our attention to the second solution (24). Calculating the necessary 
components of the superpotential, we get 

W 0a = ^-(l-e-^ 13 ). (37) 

Substituting from (37) into (34) we get 

E(r) = 2m (l - e" r3/rl3 ) . (38) 

That is twice (36): it is clear from (36) and (38) that if r — > 0, E(r) — > and as r — > oo, 
E(r) — > m for (36) and E(r) — > 2m for (38). Note that E(r) > for all r: < r < oo 



5. Main results and Discussion 



In this paper we have studied the nonsingular black hole spherically symmetric solution 
in the HNM tetrad theory of gravity. The axial vector part of the torsion, a M for these 
solutions is identically vanishing. 

Two different exact analytic solutions of the HNM field equations are obtained for the 
case of spherical symmetry. The two solutions give rise to the same Riemannian metric ( 
spherically symmetric nonsingular black hole metric). The exact solutions (18) and (24) 
represent a black hole which contain the de Sitter world instead of a singularity. The stress- 
energy momentum tensor (20) responsible for geometry describes a smooth transition from 
the standard state at infinity to isotropic state at r — > through anisotropic state in inter- 
mediate region. This agrees with the Poisson-Israel prediction concerning "non-inflationary 



material at the interface" JO] 



It was shown by M0ller |3] that a tetrad description of a gravitational field equation 
allows a more satisfactory treatment of the energy-momentum complex than does general 



relativity. Therefore, we have applied the superpotential method given by Mikhail et al.flo 
to calculate the energy of the central gravitating body. It is shown that the two solutions 
give two different values of the energy content. The following suggestions may be considered 
to get out of this inconsistency: 

(a) The energy-momentum complex suggested by M0ller 0] is not quite adequate; though 
it has the most satisfactory properties. 

(b) Many authors believe that a tetrad theory should describe more than a pure gravitation 
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field. In fact; M0ller himself || considered this possibility in his earlier trials to modify 
general relativity. In these theories, the most successful candidates for the description of 
the other physical phenomenon are the skew-symmetric tensors of the tetrad space, e.g., 
— The most striking remark here is that: All the skew-symmetric tensors vanish 



for the first solution; but not all of them do so for the second one. Some authors; e.g; p2[ |3 



believe that these tensors are related to the presence of an electromagnetic field. Others; 



e.g.; [£4|] believe that these tensors are closely connected to the spin phenomenon. There 
are a lot of difficulties to claim that HNM theory deserves such a wider interpretation. This 
needs a lot of investigations before arriving at a concrete conclusion. 

(c) Other possibility is that HNM theory is in need to be generalized rather than to be rein- 
terpreted. There are already some generalizations of HNM theory. M0ller himself considered 
this possibility at the end of his paper ||; by including terms in the Lagrangian other than 
the simple quadratic terms. Sdez |7j] has generalized HNM theory in a very elegant and natu- 
ral way into a scalar tetradic theories of gravitation. In these theories the question is: Do the 
field equations fix the tetradic geometry in the case of spherical symmetry? This question 
was discussed in length by Sdez p5| . The results of the present paper can be considered as 



a first step to get a satisfactory answer to this question. Mayer H has shown that M0ller's 



theory is a special case of the poincare gauge theory constructed by Hehl et al. pOfl . Thus 
poincare gauge theory can be considered as another satisfactory generalization of M0ller's 
theory. 



(d) Mikhail et al. [15] calculated the energy of two spherically symmetric solutions and 
found that the energy in one of the two solutions does not coincide with the gravitational 
mass. Shirafuji et al. ]16| extended the calculations to all the stationary asymptotically 



flat solutions with spherical symmetry, dividing them into two classes, the one in which the 
components (Ao) and (A«) of the parallel vector fields (Am) tend to zero faster than l/\/r 

for large r and the other, in which those components go to zero as 1/ \fr. It was found that 
the equality of the energy and the gravitational mass holds only in the first class. It is of 
interest to note that the two tetrad structures (21) and (25) have those properties, i.e., the 
first tetrad structure (21) the components (Ao) and (A a ) go to zero faster than for 



large r and so its energy is the same as that obtained before [p7 |. As for the second tetrad 



structure (25) the components (Ao) and (A tt ) go to zero as l/VF. So its energy content is 



different from the energy content of the first solution (18) and from that given by Yang |26 
and Radinschi |27| by factor 2. 



It is of interest to note that we have obtained two exact solutions for the partial dif- 
ferential equations (16) under some special constraint. The general solution for the partial 
differential equation in the case when the energy-momentum tensor does not vanish is not 
yet obtained. This will be studied in future work. 
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